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Abstract 

The electromagnetic elastic form factors of pseudoscalar and vector mesons are ana- 
lyzed for space-like momentum transfers in terms of relativistic quark models based on 
the Hamiltonian light-front formalism elaborated in different reference frames (q + = 
and q + ^ 0). As far as the one-body approximation for the electromagnetic current 
operator is concerned, it is shown that the predictions of the light-front approach at 
q + = should be preferred, particularly in case of light hadrons, because of: i) the 
relevant role played by the Z-graph at q + ^ 0, and ii) the appropriate elimination of 
spurious effects, related to the orientation of the null hyperplane where the light-front 
wave function is defined. 



PACS numbers: 12.39.Ki; 14.40. Cs; 13.40. Gp 

Keywords: Electromagnetic form factors; Relativistic quark model; Light-front 
formalism 



To appear in Physical Review C. 



1 Introduction 



The Hamiltonian light-front (LF) formalism is one of the most popular techniques developed 
to treat relativistic bound states of systems containing a fixed number of constituents [|IJ. The 
LF formalism is characterized by the fact that it maximizes the dimension of the subgroup 
of kinematical (interaction free) generators of the Poincare group. As a matter of facts, 
the interaction term appears only in three out of the ten Poincare generators, namely in 
the "minus" component of the four-momentum (P~ = P° — P z ) and in the two transverse 
rotations about the x and y axes, with the null-plane being defined by x + = t + z = 0. 

Hamiltonian LF quantum models have been widely used for the investigation of 
hadronic form factors, like, e.g., the electroweak form factors of mesons and baryons within 
the framework of the constituent quark picture of hadrons, or the elastic electromagnetic 
(e.m.) form factors of the deuteron viewed as a composite two-nucleon system. A relevant 
issue is to know to what extent the response of a composite system can be understood in 
terms of the properties of its constituents. To this end the one-body approximation for the 
current operator has been extensively considered. In case of spin- 1/2 constituents and for 
the e.m. current operator, which are of interest in this work, one has 
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where q 2 = q ■ q is the squared four-momentum transferred to the system and f^L is the 
Dirac(Pauli) form factors of the j-th constituent with mass rrij. In this paper we will limit 
ourselves to space-like q (i.e., q 2 < 0). 

While the full current J M is covariant with respect to the (interaction dependent) 
transverse rotations, its one-body approximation ([!]) is not, and the most direct manifestation 
of the loss of the rotational covariance is the so-called angular condition. Indeed, it is well 
known |0J that all the form factors appearing in the covariant decomposition of a conserved 
current can be expressed in terms of the matrix elements of only one component of the 
current, namely the plus component J + = J° + J z . It may occur however that the number 
of form factors is less than the number of independent matrix elements of the plus component 
obtained from general properties of the current operator. This means that in such situations 
a relation among the matrix elements (the angular condition) should occur in order to 
constrain further their number. The use of the one-body current ([]]) may lead to important 
violations of the angular condition, which do not allow to extract uniquely the form factors 
from the matrix elements of jt^ (cf. Refs. |2|, |3|] for the case of the p- meson). 

Within the LF formalism two different approaches [§], || have been proposed to over- 
come the angular condition problem. Both approaches really solve the angular condition 
problem and make use of the plus and transverse components of the one-body current ([J). 
However, the approach of Ref. 0] is realized at q + = and does not introduce explicitly 
any covariant current, whereas the approach of Ref. || is developed at q + ^ and a covari- 
ant approximation of the current J M is explicitly constructed. The two approaches, which 
we stress are both based on the one-body approximation ([[]), are inequivalent, because the 
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Lorentz transformation connecting a frame where q + = to a frame where q + 7^ is inter- 
action dependent. In other words, the impact of the additional many-body currents needed 
to construct the full current J M , might be substantially different at q + = and at q + 7^ 0. 

The aim of this paper is to address the issue of the relevance of the many-body 
currents by comparing the predictions of the approaches of Refs. ||| and || in case of the 
(space-like) e.m. elastic form factors of both light and heavy pseudoscalar and vector mesons 
adopting the general framework of the constituent quark model. It will be shown that the 
two above-mentioned LF approaches are inequivalent because of the different contribution 
of the so-called Z-graph |J at q + = and at q + 7^ 0. While at q + = it is possible to cancel 
out exactly the Z-graph (see Ref. ||), the latter is active at q + 7^ 0, but ignored in Ref. [|J. 
Moreover, it will be shown that the Z-graph provides an important contribution in case of 
light hadrons, whereas it vanishes in the heavy-quark limit, where the two LF approaches 
predict the same universal Isgur-Wise (IW) function 0, ||. It will be pointed out that within 
the approach of Ref. || the form factors are function of Q 2 /M 2 , where M is the mass of 
the hadron, and such a dependence is not efficient for describing the phenomenology of light 
hadrons. Furthermore, in case of vector mesons, the spurious effects related to the orientation 
of the null hyperplane where the LF wave function is defined, can be properly eliminated 
in the LF approach at q + = ||, while they are ignored and cannot be eliminated within 
the LF approach at q + 7^ 0. Thus, the latter approach appears to maximize the impact 
of the additional many-body currents needed to achieve consistency with experiment, and 
therefore, as far as the one-body current ([I]) is concerned, the predictions of the LF approach 
at q + = should be preferred, particularly in case of light hadrons. 

The plan of the paper is simply as follows. Section 2 is devoted to a brief description 
of the two LF approaches, which are then applied to the evaluation of the elastic form factor 
of both light and heavy pseudoscalar mesons. The case of vector mesons, where the angular 
condition becomes manifest, is illustrated in Section 3. The conclusions are summarized in 
Section 4. 



2 Pseudoscalar mesons 

In this Section we consider a pseudoscalar (PS) meson with mass Mps, made of two con- 
stituent quarks q± and q~2 with mass mi and 777,2, and with electric charges e\ and e%. The 
matrix elements of the e.m. current for the elastic channel are given by 

(P'\r\P) = F PS {Q 2 ) {P + P'Y , (2) 

where Q 2 = — q 2 = —(P' — P) 2 is the squared four-momentum transfer (Q 2 > 0) and Fp$(Q 2 ) 
is the elastic form factor. For both the approaches of Refs. || and we use the same LF 
wave function, which as well known jl] can be factorized into the product of a center-of-mass 
and an intrinsic parts. In terms of the intrinsic LF variables, defined as 

e = P t/p + = i-pt/p + , 

k x = Pi±-£P± = -P2± + (1-0P±, (3) 
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where \P±, P + ) describes the LF center-of-mass state [normalized according to (P'±, P' + \P±, 
P+) = 2P+ 5(P' ± - P ± ) 5(P' + - P+)}, A(£, k±) = M Q [1- (mf - m 2 2 ) 2 /M 4 ]/4£(l - f ) is a 
normalization factor for the intrinsic LF wave function, wps{k) is the radial wave function 

(normalized as / °° dk k 2 w 2 PS (k) = 1) with k = yjk 2 ± + k 2 , k z = M (£-l/2) + (m 2 -m 2 )/2M , 
and Mq is the free mass given explicitly by 



M, 



ml + kx , Tn% + k 2 



+ 



1-e 



(5) 



Finally, the quantity R^ PS ^ appearing in Eq. (|J) is the product of (generalized) Melosh 
rotation spin matrices |ll]|, viz. 



R {PS \t,h 



(A 2 |i4(l-£,-£ ± ,m 2 )|A' 2 > 



with 



RM^,k±,m 1 ] 



mi + ^Mo — «<5 • (z x 
(mj + £M ) 2 + A;i 



(6) 



(7) 



In terms of the LF spinors w(pi, Ai) and v(j)2, A 2 ) the Melosh factor R^ PS ^ can be conveniently 
written as (cf. Refs. @(b) and §) 



Ai A2 



=u(pi, Ai)7 5 w(p 2 ,A 2 



mi - m 2 J 



As for wps(k), in what follows we will consider a specific choice, namely the eigen- 
functions of the quark potential model of Ref. because the latter nicely reproduces 

the mass spectra of both light and heavy mesons, which are of interest in this workf]. The 
constituent quark masses used throughout this paper are taken from Ref. fll2f , namely: 
m u = md = 0.220 GeV, m s = 0.419 GeV, m c = 1.628 GeV and m b = 4.977 GeV. 
The masses of the corresponding charged PS and vector mesons are taken from PDG [|13| . 
namely: M w (M p ) = 0.1396 (0.767) GeV, M K (M K *) = 0.4937 (0.892) GeV, M D (M D *) = 
1.869 (2.010) GeV and M B (Mb*) = 5.279 (5.325) GeV. 

b Actually different choices for wps(k) are clearly possible, but the qualitative results presented in this 
paper will not change. 
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Let us now briefly describe the basic features of the approaches of Refs. |4|] and [0 . In 
the former approach the working frame is a Breit frame where q + = and the e.m. current 
operator is given by the one-body approximation (JJ). However, the matrix elements of the 
one-body current do not have the decomposition given by Eq. (0). This is related to the 
fact that the LF wave function (^) is specified on the null hyperplane, whose orientation can 
be identified by its normal four-vector u. The standard choice of uj is along the minus axis. 
Note that uj is a null four- vector (uj ■ uj — 0) and uj ■ q = q + . 

As firstly pointed out in Ref. and subsequently derived from an analysis of the 
Feynmann triangle diagram in Ref. [f§]f| the matrix elements of an approximate current do 
depend upon the four- vector uj, while only the amplitudes of the full current are completely 
independent of uj [see Eq. @]. Then the general structure of the matrix elements of the 
one-body current ([!]) for an elastic process involving a PS meson is given by [0, |J 



lf{ p V » \p )lf = f ^) ( P + py + b $Uq 2 ) —= , (9) 

UJ ■ T 

where Bpl(Q 2 ) represents a spurious form factor. Since the standard choice of u is along 
the minus axis, the physical form factor Fpg(Q 2 ) can be obtained using the plus component 
of the one-body current ([!]). Adopting the usual Breit frame where q + = q~~ = 0, one has 



P'+ = P+ = \JM PS + Q 2 /4 and P' ± = -P ± = q ± /2 with Q 2 = q\. Assuming for sake of 
simplicity point-like constituents, the elastic form factor Fpg(Q 2 ) can be cast in the following 



form 



where @, |, 10 



F«(Q 2 ) = e 1 H 1 {Q 2 ) + e 2 H 2 (Q 2 ) , (10) 

wps(k) wps(k') 

(11) 



H^Q 2 ) = j^dij dk ± \/A(C,k ± ) A(£,k'_ 



lo J " " " 1 ;•; 



V 2 (0 + kl^ 2 (0 + k'{ 

with /i(^) = mi (1 — + m 2 ^ and k'± = k± + (1 — £)q±. The explicit expression for the 
form factor Fl2(Q 2 ), corresponding to the coupling of the virtual photon with the antiquark 
q~2, can be easily obtained from Eq. ([□]) by using simply k'j_ — k± — £q±. Note that the 
form factor ( [TOD is a function of Q 2 and of the constituent masses mi and m^. 

In case of the approach of Ref. the working frame is a special Breit frame where 
the four-momentum transfer q is along the spin-quantization axis, i.e. the z axis. Therefore 

c The appearance of the four-vector 10 in the Feynmann triangle diagram can be understood as due to 
the need of expressing the off-mass-shell constituent momenta, which naturally appears in any Feynmann 
diagram, in terms of the on-mass-shell constituent momenta, which characterize the LF formalism (see Ref. 

§)• 
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one has u ■ q = q + 7^ 0, more precisely: q + = — g = Q, q± = 0, P + = J M PS + Q 2 /4 — Q/2 



P' + = y Mp S + Q 2 /4 + Q/2 and P± = P'± = 0. In such a special frame and in case 
of elastic processes a choice for the e.m. current operator compatible with (extended) 
Poincare covariance and hermiticity is given by 

f = \{ C " + V"[rx{-*)] e lwSx C* u e^ 3 '} , (12) 

where r x {— 7r) represents a (— 7r) rotation around the x axis, is the x-component of the LF 
spin operator and L(£) is the element of the Lorentz group corresponding to £ e SL(2, C). In 
Eq. (12) the operator C M should fulfill the (extended) Poincare covariance and the specific 
choice made in Ref. || is as follows: 



C+ = 


J+ 

J (l) 


c± = 


J(l) 


c~ = 


J+ 
J (l) 



(13) 

where the last equation ensures the gauge invariance of the current operator j M . The covariant 
decomposition of LF(P'\j fJ '\P)LF is therefore given by 

LF (P>\f\P)LF = F p l) S {P + Pr , (14) 

where no spurious structure is present thanks to the gauge invariance of and to the fact 
that uj ■ q = q + ^ 0. Then, the elastic form factor Fpg can be calculated using the plus 
component j + , whose matrix element LF(P'\j + \P) lf reduces to the corresponding one of 
the plus component of the one-body current Jt^ evaluated at q + = — q~ = Q (see Ref. ||). 
Therefore, one gets 

= e 1 H 1 (Q 2 /M 2 PS ) + e 2 H 2 (Q 2 /M 2 PS ) (15) 

with 

H 1 {Q 2 /M 2 PS ) = ^— K - £ d£ J dk ± s/a(£, k ± ) A(?, k ± ) 



1 - k/2 

wps(k) w P s(k') /i(0/i(O + k 2 ± 



(16) 



where = k + (1 — k) £ and 



q + _ Q/M P s Q 4 Q 2 Q 2 

P' + Jl + Q 2 /AM 2 PS + Q/2M PS V 4M 4 5 M PS 2M PS ' 



The explicit expression for the form factor H 2 {Q 2 /M 2 PS ) can be easily obtained from Eq. 
(Il6|) by using simply £' = (1 — k) £. Note that the factor (1 — ac)/(1 — k/2) in the r.h.s. of 
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Eq. fllTf ) is nothing else than the factor 2P + /(P + + P f+ ), where the numerator comes from 
the normalization of the LF center-of-mass states. Moreover, one has < k, < 1. 

In Eq. (fl5|) we have explicitly taken into account that within the approach of Ref. 
H the elastic form factor J-'pg is not a function of Q 2 , but of the ratio Q 2 /Mp s (and of the 
constituent masses m\ and m2). This is an important feature, which naturally emerges when 
a " longitudinal" frame with q + ^ and q± = is chosen. As a matter of facts, since only the 
LF fraction £ changes in the final state, the form factor J^pg cannot be function of q + only, 
but of q + /P + (or equivalently q + / P' + = k). Thus, the dependence of the elastic form factors 
upon the ratio between the momentum transfer and the mass of the system is expected to 
characterize the predictions of the LF approach of Ref. f| for a generic hadron^ and this 
fact has important consequences. The most direct one is that the charge radius of a hadron 
is expected to be approximately proportional to the inverse of its mass. Such a behavior is 
completely at variance with the phenomenology of light hadrons. Indeed, for instance, the 
experimental charge radii of the pion, the kaon and the nucleon are approximately the same 
[rl h = 0.660 ± 0.024 fm [0(a), r* h = 0.58 ± 0.04 fm |5[(b) and r p ch = 0.883 ± 0.014 fm 
||15|| (c)1, while the kaon (nucleon) mass is larger than the pion mass by a factor of ~ 3.5 
(6.7). 

From the above considerations it is quite clear that the approach of Ref. || is likely 
to be not efficient for the description of the phenomenology of light hadrons; in other words, 
substantial effects of additional many-body currents should be explicitly considered in order 
to achieve consistency with experiment. An explicit demonstration is provided by the pion 



case. The results obtained using Eqs. ( p^|16| ), are reported in Fig. 1 and compared with 
the corresponding results of the approach of Ref. |4j] based on Eqs. ( p!0| - |l~T|) . The form factor 
evaluated at q + ^ exhibits a very rapid falloff with increasing Q 2 , corresponding to a charge 
radius of ~ 5 fm. On the contrary, the form factor obtained at q + = is much higher and 
corresponds to a charge radius of ~ 0.46 fm. Similar results hold as well in case of the 
kaon, whose charge radius turns out to be ~ 2 fm at q + ^ and ~ 0.43 fm at q + = 0. All 
the calculations have been done assuming point-like constituents and therefore only within 
the approach at q + = it is possible to recover the agreement with the experimental charge 
radius of the pion (and the kaon) by introducing a constituent size of ~ 0.45 fm, as already 
proposed in Ref. JK|. The same constituent size is suggested also by the recent analysis of 
the elastic nucleon form factors carried out within the covariant LF approach at q + = in 



Ref. 16 



The question to be addressed now is clearly the origin of the large differences in the 
pion form factor evaluated at q + = and q + ^ 0. The answer is already well known from 
the works of Ref. || and more recently from the results of Refs. 0, [17]]. There it has been 
shown that the one-body LF form factor at q + = [Eqs. (JlO-11)] matches exactly the form 
factor embedded in the Feynmann triangle diagram, evaluated using the one-body current 
(H) and adopting the appropriate bound-state vertex corresponding to the LF wave function 
(|j) [see Eq. (11) of Ref. [HI]. The matching is due to the vanishing of the contribution of 

d This means that the LF approach at q + ^ provides naturally implicit many-body contributions 
associated with the hadron mass. 
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Figure 1. Elastic form factor of the pion, F n (Q 2 ), versus Q 2 . Solid line: LF approach of Ref. 
0] at q + = 0, corresponding to Eqs. (|iO|-|ll|). Dashed line: LF approach of Ref. || at q + ^ 0, 
corresponding to Eqs. ( |l5| - [T^ ) ■ Point-like constituent quarks are assumed in the calculations. The 
radial wave function wps(k) is taken to be the eigenfunction of the quark potential model of Ref. 

B 



the Z-graph at q + = 0. Thus, one has Fp$ = p^ ian 9 le \ Thanks to the covariance property 
of the triangle diagram, the same form factor can be obtained in a frame where q + ^ 0. 
However, in such a frame the anatomy of the form factor is different in the sense that it is 
given by the sum of two non-vanishing contributions, the spectator and the Z-graph terms: 
p(tmangie) _ jri^v-) _|_ j^^~9 ra P h ) ^he spectator term (see Ref. |I7||), evaluated in the Breit 
frame where q + = —q~ = Q, turns out to coincide with Eqs. (|T5| - p!6D , i.e. with the form 
factor predicted by the approach of Ref. 0: Tp S = Fps^ ■ Thus, finally one has 

fps-Sps + fps J { L °J 

where J r P Z s ~ graph ' > is the contribution of the Z-graph in the triangle diagram evaluated in 
the Breit frame where q + = —q~ = Q. Thus, the origin of the differences in the pion form 
factor evaluated within the LF approaches of Refs. Q and || is the Z-graph contribution at 
q + 7^ 0, which is ignored in Ref. |J. From Fig. 1 it is clear that the Z-graph dominates the 
pion form factor at q + ^ 0, while the contribution of the spectator term is almost negligible 
(except at the photon point). Therefore, since the Z-graph is a many-body process, we can 
conclude that the approach of Ref. ||, based on the specific choice (|T3| ) and on the Breit 
frame where q + = —q~ = Q, appears to maximize the impact of the many-body currents 
needed for consistency with experiment, particularly in case of light hadrons. 



S 



Since the Z-graph is expected to vanish in the heavy-quark limit, it is worthwhile to 
study the behavior of the two LF approaches in that limit. We expect that for fixed values 
of the dot product w of the initial and final meson four-velocities, given by 

w = P- P'/M 2 PS = 1 + Q 2 /2M 2 PS , (19) 

one should have 

linVni^oo 

Hi(Q 2 ) = lim 

mi— >oo Hx(Q 2 /M 2 PS ) = Ciw(w) , (20) 

where £iw( w ) is the Isgur-Wise (IW) form factor J7|, which has been already calculated 
within the LF approach at q + = in Ref. [|18|]. We have therefore calculated the form 
factors H^Q 2 ) [Eq. flTTD] and Hi(Q 2 /Mp s ) [Eq. flTE|)] for various values of the constituent 
mass mi pertaining to the cases of ir, K, D and B mesons, keeping fixed the constituent 
mass m2 at the value 7712 = 0.220 GeV. The results are reported in Fig. 2 in terms of the 
variable w and compared with the IW function ^iwiw) calculated in Ref. jl 




w 



w 



Figure 2. The PS form factors H\ [Eq. (|IT|)1 (a) and Ji\ [Eq. (0)] (b) versus the variable w 
[Eq. (19)]. Triple-dot-dashed, dotted, dot-dashed and dashed lines correspond to the cases of n, 
K, D and B mesons, respectively. In (a) and (b) the solid line is the same IW function (,iw(w) 
as calculated in Ref. |18|| . In (b) the dashed line almost coincides with the solid line. Point-like 
constituent quarks are assumed in the calculations. The radial wave functions wps(k) are taken to 
be the eigenfunctions of the quark potential model of Ref. [12]. 



Few comments are in order: i) Eq. ( pO"D is fulfilled, i.e. the two LF approaches 
predict the same asymptotic IW function. This confirms that the differences between the 
results at q + = and q + 7^ are entirely due to the Z-graph contribution ignored at q + 7^ 0; 
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ii) the convergence of the calculated form factors to the asymptotic IW function occurs from 
below at q + 7^ and from above at q + = and it is much faster in the former case: TCx ~ 
already at the B-meson mass (see Fig. 2(b)); iii) the dependence of the calculated form 
factor at q + 7^ upon the constituent mass mi is quite mild at q + 7^ 0, which implies that in 
terms of the variable w the approach of Ref. H predicts quite similar form factors for light 
and heavy mesons. Such a prediction is not reasonable; indeed, the dynamics of a heavy 
quark is characterized by the Heavy Quark Symmetry (HQS) [0, |§, which makes a heavy 
quark blind to the spin and flavor of light spectator quarks. The situation is opposite in light 
hadrons, where already the mass spectrum clearly exhibits a rich spin and flavor dependent 
structure. Thus, at q + 7^ the form factors of light hadrons are basically dominated the 
Z-graph, while the one-body spectator term [i.e., Eqs. p!5| - |l6D l plays only a marginal role. 
On the contrary, at q + = the Z-graph is suppressed and the main contribution to the form 



factors is given by the one-body spectator term [i.e., Eqs. 10-11)] . Therefore, the approach 
of Ref. || should be improved by making a choice of the operator C M different from Eq. 
(fL3|). If this will be done in such a way to include the effects of the Z-graph, we expect that 
the present one-body predictions at q + = will be recovered as the sum of one-body and 
additional many-body contributions at q + 7^ 0. 

If we scale the result obtained for r^ h at q + 7^ (i.e. r* h ~ 5 fm) according to the \ jM 
dependence, where M is the mass of the system, one may obtain a charge radius of ~ 0.7 fm 
for M ~ 1 GeV, i.e. around the nucleon mass. In other words it might happen that the 
approach of Ref. could reproduce the proton charge radius even assuming point-like (or 
almost point-like) constituents. Moreover, thanks to the approximate dipole behavior of the 
proton form factors at least for Q 2 < M 2 ~ 1 GeV 2 , a good agreement with the experimental 
data at low Q 2 might be achieved at q + 7^ without the need of introducing a finite size for 
the constituent quarks. However, it should be clear that such an agreement is not physically 
meaningful unless one demonstrates that the Z-graph gives a negligible contribution. This 
is likely not to be the shown in Fig. 3, where the quantity Rz, defined as 

Rz = \ , (21) 
All 

is reported as a function of the variable w for the various pseudoscalar mesons so far consid- 
ered. The quantity Rz is a measure of the relative importance of the Z-graph with respect 
to the calculated form factor at q + 7^ 0. Since Mk is almost half of the nucleon mass and 
Mb approximately twice the nucleon mass, we expect that, as a conservative estimate, the 
Z-graph contribution for a system with a mass around 1 GeV should be between the dotted 
and dashed lines shown in Fig. 3, and therefore not negligible. 

The above criticisms can be directly extended to the recent results obtained in Ref. 
19fl for the nucleon elastic form factors within the so-called Point Form Spectator Approx- 



imation (PFSA). Such an approach was firstly proposed and illustrated in Ref. [^UJ; here 
it suffices to say that it is based on the one-body approximation ([I]) for the e.m. current 
operator, but carried out within a form of the dynamics different from the LF one, namely 
the point form. However, as in the case of the LF approach at q + 7^ 0, within the PFSA the 
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Figure 3. The ratio Rz [Eq. (|2l|)l obtained using the results shown in Fig. 2, versus the variable 
w [Eq. fll9|)1. The meaning of the lines is the same as in Fig. 2. The vertical long-dashed line at 



w 



1.5 corresponds to Q 2 = M 2 , where M is the mass of the system. 



elastic form factors of a hadron are not function of Q 2 , but of Q 2 /M 2 , where M is the mass 
of the hadron (cf. Ref. [p0[|). Moreover, within the PFSA the Z-graph is totally ignored 
(i.e., no match with the Feynmann triangle diagram). Therefore, the agreement with the 
elastic nucleon data obtained in Ref. |L9| assuming point-like constituent quarks, is likely to 



be physically meaningless in the same way and with the same arguments already explained 
in case of the LF approach at q + ^ 0. A more detailed analysis of the PFSA is in progress. 

Before closing this Section, we want to stress again the main result so far achieved, 
namely: among various Hamiltonian formalisms (point as well as light-front forms) based 
on the one-body e.m. current ([[]), the predictions of the LF approach at q + = should be 
preferred, particularly in case of light hadrons. This is due to the fact that when q + ^ the 
relevance of many-body currents appears to be amplified by the occurrence of the Z-graph 
term, which turns out to be essential also to compensate the unwanted Q 2 /M 2 dependence 
of the one-body form factors calculated in " longitudinal" frames. Finally, we point out that 
the choice q + = is possible only for space-like q (q 2 < 0). Indeed, for time-like q (q 2 > 0) 
one has always q + ^ 0. In this case one needs to perform an analytic continuation from 
space-like to time-like q. This cannot be easily done in the standard LF formalism because 
the contribution of the Z-graph cannot be eliminated when q 2 > (see Refs. \I7[ 



However, a proper analytic continuation can be achieved for the Feynmann triangle diagram 
by means of the so-called dispersion approach, which is described in Ref. [|^] and has been 



extensively applied to time-like processes, like heavy meson weak decays, in Ref. [23]. We 
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stress that for space-like q the dispersion approach result matches the LF one at q + = (see 
Ref. H). 



3 Vector mesons 

The elastic e.m. response of a vector system is described by three physical form factors 
Fi{Q 2 ) {i = 1)2,3), which appear in the following covariant decomposition of the matrix 
elements of the e.m. current operator J M : 

JV, 8 ) = -{P + P'T U(Q*)e*(P>, s') ■ e(P, a) + [e*(P', s>) ■ q] [e(P, s) ■ g] J 

+ F 3 (Q 2 ) {[e»(P', s'))* [e(P, s) ■ q) - e»(P, s) [e*(P', s') ■ q)} , (22) 

where e(P, s) is the LF polarization four-vector of the spin-1 system (with mass My) corre- 
sponding to spin projection s and total four-momentum P. 

In the standard LF approach at q + = (see Ref. @) only the matrix elements of 
the plus component of the one-body approximation ([I]), which will be denoted by J^(s', s), 
are taken into account. After considering general properties of the current operator (like, 
e.g., the time reversal symmetry) the number of independent matrix elements J^(s',s) 
turns out to be four, while the physical form factors are three. A further condition arises 
from the rotational invariance of the charge density, which however involves transformations 
based upon Poincare generators depending on the interaction. Such an additional constraint, 
known as the angular condition |24j], reads as 

(1 + 277) J+(l, 1) + J+ (1, -1) - v^Jjj(l, 0) - J+(0, 0) = , (23) 

where 77 = Q 2 /4M V . The angular condition fl23|) is not satisfied by the matrix elements 
J(rj(s', s ) an d therefore the extraction of the one-body form factors F^ l \Q 2 ) is not unique. 

In Ref. [|J the problem of the violation of the angular condition (^) is solved by 
considering the matrix elements of the tensor T/^ , which are related to those of the one- 
body current by 

Jfa(s l ,s) = e' a (P',s l )TSf e p (P,a). (24) 
In terms of spin-1 LF wave function one has 

T$f= LF (P'^\J^\P,(3) LF , (25) 
where for an 5-wave vector system 



P,P) LF = R( y \£,k ± ;P) w v (k) 



Atx 



(26) 
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with 



[R^(£,^;/3) 



A1A2 



\/2 ^Mq 2 - (mi - m 2 ) 2 
u(pi, Ai 



p _ (gi - ^ 
M + m 1 + m 2 



(27) 



Following Refs. |3J and |3|, 13] the general decomposition of the tensor T/|\ reads as 



(2? 



where the tensor J/^ 3 is independent of the four-vector cu, while B^?^{u) contains all the 
possible covariant structures depending on ui. One gets || £§ 



T fj,,a3 



+ 



pap/3 piap/f3 p/apfj p p 



M 2 



M 2 



+ 



2M V 



9 



pllipia' 



Ml 



pi 



- 9 



»0 



pvpP\ 



M 2 ) 



P' 



P'-q 



My 



pmpp ^(l)(g2^ 



My 



+ 



2M y 



M 2 My 



(fP? - q^P" 



+ ^ 1} (g 2 ) (g^P 13 + ^P' a ) + Hf\Q 2 )q 



(29) 



and 



2(u> • P) 



„a, .3 „3, ,a 

+ B {1) (0 2 ) q q 
+ ^ 4 {Q } 2(u ■ P) 

,8 „8, ,a 

+ B^ qu qu 



+ {p + p 



r 2 nfl)/^ <*> ^ 



M 2 ^^(g 2 ) 



2(w • P) 



J^B?(Q 2 ) 



(u-P) 



sT- , 



,pW ( o 2) 

+ ^ 8 w j q 2(u ■ P) 



(30) 



where all the covariant structures included in Eqs. (|29| - p0|) satisfy both parity and time 
reversal symmetries. 
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In Eq. ( p9|) there are seven form factors, namely the three form factors F^ 1 \Q 2 ) 
(i = 1,2,3) and the four form factors Hf\Q 2 ) (j = 1,...,4). The form factors F^\Q 2 ) 
appear in covariant structures which are transverse to all the external momenta P, P' and q, 
while the form factors H^IQ 2 ) describe the loss of transversity (including the possible loss 
of gauge invariance) of the tensor (|29|) . Therefore, since e(P, s) ■ P = e(P f , s') ■ P' = 0, the 
form factors H^\Q 2 ) do not appear in the decomposition of the matrix elements J^(s', s). 

In Eq. © all the B^ ] (Q 2 ) (k = 1,2, ...8) are spurious form factors, which can 
contribute to the matrix elements J^(s,s'), namely: 

JfaW, *) = ~(P + P'T ^Fi X \Q 2 )e\P', s') • e(P, s) + [e*(P', s>) ■ q] [e(P, s) ■ q] 

+ P 3 (1) (Q 2 ) {[e»(P', s')]* [e(P, s)-q)- e"(P, s) [e*(P', s') ■ q)} 
+ e* a {P',s')B^ e p (P,s). (31) 
For the angular condition one has || 

(l + 2r / )J+(l,l) + J+(l,-l) - 0^,7+ (1,0)- 4,(0,0) = 

= -B^iQ 2 ) - B?\Q 2 ) ^ , (32) 

which implies that at q + = the loss of rotational covariance of the one-body current ([!]), i.e. 
the violation of the angular condition (^), is described by the spurious structures containing 
B^\Q 2 ) and Bj 1 \Q 2 ) in Eq. Q3U|). Finally, note that the loss of gauge invariance of the 
one-body current ([!]) at uj ■ q = q + = is described by the spurious form factor Bq 1 \Q 2 ) 
appearing in Eq. (|30D . 

As explained in Ref . , in the Breit frame where q + = the physical one-body form 
factors F t {1 \Q 2 ) (i = 1,2,3) can be obtained through the following equations 



rn-rm 

i (l) 
IP 



2P+ ' 

rp+,yy _ rp+,XX rp+, + + 



2 V 2P+ 2(1 + v ) 2P+ y^IT^) 2P+ 



FPm = Zk. _ ^±IZk_ . ,33) 

At variance with the PS case the components of the tensor (^) with fi = y are essential for 
the extraction of the physical form factors, more precisely for the determination of F^\Q 2 ). 
Note also that the form factors F± (Q 2 ) and F^ (Q 2 ) are determined only by the components 
with /i = +. 

We stress that by means of Eq. ([!3|) the extraction of the physical form factors 
F < > l \Q 2 ) is not plagued at all by spurious effects, including those related to the loss of 
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rotational covariance and gauge invariance. The angular condition problem (^) is therefore 
completely overcome without introducing explicitly any covariant current operator. 

As far the approach of Ref. || is concerned, the matrix elements j^{s' ', s) of the covari- 
ant current (|T2|) can be expressed in terms of the LF wave function (|26"l - f27| ) as: j fM (s', s) = 
(P', s'|j M |P, s) with \P,s) = ep(P,s) \P,/3)lf- In Ref. || it is assumed that the matrix 
elements j M (s', s) have the following decomposition 



P + P'Y {J r i 1) e*{P',s')-e{P,s) + 



(i> 



2Ml 



[e*(P',s')-q] [e(P,s)-q] 



+ Jf } {[e»(P', s')}* [e(P, s).q]- e?(P, s) [e*(P', s') ■ q}} 



(34) 



and it is shown that there are only three independent matrix elements j ,x (s , ,s), namely 
ioo; 3u an d JioB- Thanks to the choice ([13]), the latter are very simply related to the matrix 
elements of the one-body current (JX|) evaluated at q + ^ 0, which will be denoted by J^{s' , s) 
to distinguish them from the matrix elements J^Js', s) evaluated at q + = 0. Ignoring for the 
moment the possible presence of spurious structures, in the Breit frame where q + = —q~ = Q 
one has 



Joo 

i x 
J 10 



^(t)(H) 



2M V 
2M V 



1 + 77 {(1 + 2rj) - 2rj (1 + rj) J- 2 (1) - 2 V ^ 3 (1) } 



1 + 77 F{ 



(i) 



^i)(10)-J ( i)(01) 



2Mv 



1 + 7] Jr)/2F. 



(i) 



(35) 



Thus, the form factors 
relations 



1,2,3) may be uniquely determined through the following 



^\Q 2 /m 2 v ) 



n l \Q 2 /Mi 



+ 



2MvVl + V (1) 




1 { 


4M V 7] 


(i + v) 1 




^(i)(10)- 




1 



2M vy /TTvV^V 



^i)(io)- J5)(oi) 



(36) 



where we have explicitly taken into account that the form factors J-[ are not function of 
Q 2 , but of the ratio Q 2 /My, thanks also to the fact that the matrix elements J~^(s',s) are 
proportional to 2P + because of the normalization of the LF center-of-mass states. Note 
that the form factor J 7 ^ is determined through the x component of the one-body current, 
while J~i depends on the plus component only. 

e Note that the matrix element jf is not independent from jf because one has j\ { 
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1 3w 



i. 
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The predictions of the LF approaches at q + = [Eq. fl53|) l and at q + ^ [Eq. (|36|)1, 
obtained in case of the p meson adopting for w p (k) the corresponding eigenfunction of the 
quark potential model of Ref. |l2 |, are compared in Figs. 4-6 in terms of the conventional 
charge Gq(Q 2 ), magnetic G\{Q 2 ) and quadrupole (j2(<3 2 ) form factors, defined as 



G (Q 2 ) 
Gi(g 2 ) 

G 2 (Q 2 ) 



Fs(Q 2 ) , 



Fs(Q 2 ) 



v)F 2 (Q 2 



^i(Q 2 )-^(Q 2 )-(i + r/)F 2 (Q 2 



(37) 



o 



1 .0 
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Figure 4. Charge form factor Gq(Q 2 ) for low (a) and high (b) values of Q 2 . The solid line is 
the result obtained within the LF approach at q + = [Eq. (|33|)1, while the dashed line is the 
prediction of the LF approach at q + ^ [Eq. (^)] . Point-like constituent quarks are assumed in 
the calculations. The radial wave function w p (k) is taken to be the corresponding eigenfunction of 



the quark potential model of Ref. \12\. 



It can be seen that: i) the charge radius, r ch = y — 6dG (Q 2 ) / (1Q 2 \q2 =0 , turns out 
to be r c h = 1.1 fm at q + ^ and r c h = 0.57 fm at q + = 00; ii) the charge form factor 
calculated at q + = exhibits a node at Q 2 ~ 3.4 (GeV/c) 2 at variance with the result 
obtained at q + ^ 0; iii) the magnetic moment, fi p = Gi(Q 2 = 0), is remarkably larger at 
q + 7^ [lip = 11.6] than the corresponding one at q + = [ji p = 2.35], which is quite close to 

■^Taking into account a constituent size of ~ 0.45 fm as in the pion and kaon cases, the LF approach at 
q + = predicts a charge radius of the p meson equal to ~ 0.75 fm. 



16 




the non-relativistic limit \x p = 2 (cf. Ref. Q); moreover, the magnetic form factor Gi(Q 2 ) 
calculated at q + ^ has a sharp upturn near the photon point completely at variance with 
the result obtained at q + = 0; iv) the quadrupole form factor G^Q 2 ) is predicted to be 
quite small at q + ^ both at low and high values of Q 2 . All these findings are a direct 
manifestation of the lack of the Z-graph in the LF approach at q + ^ 0. Thus, for the p 
meson the Z-graph plays an important role both at low and high values of Q 2 . We want to 
point out that the node obtained in the charge form factor Go(Q 2 ) at q + = may be (at 
least partially) related to the nature of the spin-spin term of the quark potential model of 



Ref. ||12|| , which is repulsive in the triplet spin state at short interquark distances. 

Let us now consider the heavy quark limit in which only one quark is coupled to the 
virtual photon and its mass goes to infinity, i.e. mi — » oo. Thanks to the HQS we expect 
that for fixed values of the variable w = P ■ P' /My = 1 + Q 2 /2M V one should have 

lim^^oo G (Q 2 ) = £iw(w) , 
lim mi ^oo Gi(Q 2 ) = ^iw{w) , 

lim^-xx, G 2 (Q 2 ) = 0, (38) 

where £iw( w ) is the same IW function encountered in the PS case. We have calculated 
the form factors G{ (i = 1,2,3) for various values of the constituent mass m\ pertaining to 
the cases of p, K*, D* and B* mesons, keeping fixed the constituent mass m-i at the value 
m2 = 0.220 GeV. The results are reported in Figs. 7-9 in terms of the variable w and 



compared with the IW function £,iw( w ) calculated in Ref. fT8|j . 

It can clearly be seen that Eq. ( p^) is fulfilled, confirming in this way that also in 
case of vector mesons the differences between the results at q + = and q + ^ are due to 
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the lack of the Z-graph contribution at q + ^ 0. We want to stress that the universality 
of the calculated IW function, implied by the fulfillment of Eqs. fl2T| ) and (|35[), is in nice 
agreement with the predictions of the HQS, which is an exact symmetry of QCD in the 
heavy-quark limit. Finally, from Figs. 8-9 a couple of comments are in order: i) the 
difference in the calculated magnetic moment is still remarkably sizable at the D* mass (i.e., 
around My ~ 2 GeV), and ii) the lack of the Z-graph at q + ^ makes the quadrupole 
form factor G 2 very small even at quite large values of the (active) quark mass. 

Let us now consider the question whether the decomposition (|34D is complete without 
introducing spurious structures analogous to those appearing at q + = in Eqs. ( |5D| - |5TD . It 
is clear that the dependence of e.m. amplitudes upon the four-vector ou, which we remind 
identifies the orientation of the null hyperplane where the LF wave function is defined, can 
occur independently of the value of q + . Therefore, at q + ^ one expects that the following 
general decomposition holds for the matrix elements J7(i)(s', s) of the one-body current: 



(i) 



J (1) (s>, s) = -(P + PT |^i (1) e*(P', s>) ■ e(P, s) + ^ [e*(P f , s') ■ q] [e(P, s) ■ q] 
+ {[e»(P', s')T [e(P, s)-q\- e^(P, s) [e*(P', s') ■ q)} 

5 (1) 



+ eUP',s')B?^ e p (P,s) (39) 



with 



M 2 

(1) 1 ; 2(u-P) 



+ B? q -\ + MyB? 

My V 3 (UJ-P) 2 
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w w 

Figure 7. Charge form factor Go(w) versus the variable w, calculated within the LF approach 
of Ref. Q at q + = (a), corresponding to Eq. (|33|), and the one of Ref. || at q + / (b), 
corresponding to Eq. (|36|) . Triple-dot-dashed, dotted, dot-dashed and dashed lines correspond to 
the cases of p, K*, D* and B* mesons, respectively. In (a) and (b) the solid line is the same IW 
function £,iw{w) as calculated in Ref. In (a) the dot-dashed and dashed lines almost coincide. 
Point-like constituent quarks are assumed in the calculations. The radial wave functions wy{k) are 
taken to be the eigenfunctions of the quark potential model of Ref. JO] . 



+ B 



+ 



~ti, .8 

(i) q ^ 



2(u-P 



+ B, 



(!) q u H 



qP(jj a 



2(cu ■ P) 



9 



qV-qP s 

q 2 



OJ ' 



+ (P + P 



• py 



+ 



(!) ^ + 

+ q 2( W • P) 



(40) 



where P = (P + P')/2. With respect to Eq. (|30|) we have introduced a different notation 
for the spurious form factors, because the latter may depend in general on uj ■ q, i.e. they 
can be different at q + = and q + ^ 0. Note that at q + ^ the loss of gauge invariance of 
the one-body current is described not only by the spurious structure containing B^ (as in 
the case q + = 0), but also by the ones proportional to with k = 1, ...,4. 

By means of the four- vector u the choice (|13| ) for the operator C^ 1 can be conveniently 



cast into the following form: C 7 ' = J^, — uj^(q ■ J^)/(uj ■ q). Therefore, using Eqs. (|3 



40) one immediately obtains that the spurious structures containing the form factors BfP 



with k = 1,...,4 cannot appear in the decomposition of the matrix elements C^(s',s) of 
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Figure 8. The same as in Fig. 7, but for the magnetic form factor Gi(w). 



the gauge-invariant operator C M . Consequently, also for the matrix elements j M (s', s) we 
generally expect that 



i P P>y { jfV(P', s>) ■ e(P, s) + ^' ^ ' D ' ' 



e*(P',s')-q] [e(P,s)-q] 



+ {[e»(P', s'))* [e(P, s)-q)- e"(P, s) [e*(P', s>) ■ q)} 



+ e* a (P',s')Blf e p (P,8) 



(41) 



with 



BgfH = {p + p'Y 



. ,a. .3 „ol, ,3 „3 , ,a 

My 8^ -^L-„ + £#> q U ~ q U 



( U ■ py 



2(u ■ P) 



+ 



My gfl, 



(U-P) 



+ B, 



(1) 



T 



u-q) 



q^P + qP^ 
2(uj ■ P) 



(42) 



In order to prove that the decomposition (0) is not complete, it is enough to demon- 
strate that at least one of the four spurious form factors B^_ 8 is non vanishing. To this end 
let us first consider that the possible presence of spurious structures corresponds to replace 
Eq. (M) with 



Jo + o = Jj)(00) = 2Af v v /T+^{(l + 2r / ) jf> - 27/(1 + 7?) jf> - 2 V 



20 



0.05 



0.00 



-0.05 



-0.10 - 



-0.15 - 



-0.20 - 



-0.25 




0.01 



CD 



0.00 



-0.01 



-0.02 



I 1 1 1 1 1 1 1 1 1 1 - 



(b) q + * 



_. 

^ 

\ ■■ 

- \ 



j i i I i i i I i i i_ 

2 3 



1 2 3 4 1 

W W 
Figure 9. The same as in Fig. 7, but for the quadrupole form factor G2{w). 



+ 



B 



(i) 



l + V 



+ 2 v B t 



l + V 



2 V B; 



(1) 



7« 



j+ = J+(ll) = 2M vy /l+rin 

' r J5)(10)- J5)(01)]=2M V 



j 1 
J w 



2(1 + 77) 



(43) 



Therefore, the extraction of the form factors and may be plagued by spurious 
effects, while only Ti is free from spurious effects. We want to point out that the relation 
Jio = * 3vo 1S n °t modified by the possible presence of spurious structures in Eqs. fl4"H - |4"2]) ; 
in other words the number of independent matrix elements of the current operator j M is still 
three even in presence of the spurious form factors appearing in Eq. (fj2"P . 
In addition to Eq. (@3]) 

one has 



7(1) 



J0){W) + J(i) (oi) 



Mi 



V2 



B\ L 



(44) 



which allows us to calculate the spurious form factor B^p directly in terms of matrix elements 
of the one-body current. The explicit calculation of the l.h.s. of Eq. ([44]), carried out for 
the case of the p meson, is reported in Fig. 10. It can be seen that is just a small 
fraction of J 7 ^ and therefore its impact on the extraction of T^p is quite limited. However, 
what really matters is not the quantitative impact of B^\ but the conceptual fact that it is 
non vanishing, which demonstrates that the covariant decomposition ( |34"D is not complete. 
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Figure 10. Ratio of the spurious form factor calculated from Eq. (|44|), to the magnetic form 
factor J% , obtained from Eq. ( |36| ) in case of the p meson within the approach of Ref. || at q + 7= 
(dashed line). 

Note that at the photon point B^ = 0, so that the anomalous result n P = 11.6 obtained at 
q + 7^ does not depend on spurious effects. Furthermore, the other spurious form factors 
B§~\ £>g^ and B^ entering Eq. fl43D , cannot be calculated in terms of matrix elements of 
the one-body current. Thus we can conclude that the approach of Ref. not only ignores 
the spurious terms arising from the orientation of the null hyperplane, but also it is not able 
to eliminate consistently such spurious effects in the extraction of the physical form factors. 

4 Conclusions 

In this paper we have carried out a detailed comparison of the predictions of the light-front 
approaches of Refs. (3| and || in case of the (space-like) electromagnetic elastic form factors 
of both light and heavy pseudoscalar and vector mesons adopting the general framework of 
the constituent quark model. The two approaches are based on the one-body approximation 
(|I|) for the electromagnetic current operator, but they are elaborated in different Breit frames, 
namely at q + = j|] and (| + / |. It has been shown that: i) the two light-front 
approaches are inequivalent because of the different contribution of the Z-graph at q + = 
or q + ^ 0. While at q + = it is possible to cancel out exactly the Z-graph (see Ref. Q]), the 
latter is active at q + ^ 0, but ignored in Ref. 0; ii) the Z-graph provides an important 
contribution in case of light hadrons, whereas it vanishes in the heavy-quark limit, where 
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the two light-front approaches predict the same universal Isgur-Wise function. 

We have also pointed out an important feature of the approach of Ref . || , namely: 
the elastic form factors of a hadron are basically function of Q 2 /M 2 , where M is the mass of 
the hadron. We have shown that such a dependence (shared also by the point-form approach 



of Refs. [IS, p0|) is not efficient for describing the phenomenology of light hadrons. In other 



words, the light-front approach at q + ^ maximizes the impact of the many-body currents 
needed to achieve consistency with experiment. 

Finally, in case of vector mesons, the spurious effects related to the orientation of the 
null hyperplane where the LF wave function is defined, have been analyzed in details. While 
such unwanted effects are properly eliminated in the approach of Ref. [f|], they are ignored 
and cannot be eliminated within the approach of Ref. ||. Thus, we can conclude that, as 
far as the one-body current (|I|) is concerned, the predictions of the light-front approach at 
q + = should be preferred, particularly in case of light hadrons. 
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